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I.  Introduction 


1 2 

Various  methods  ’ have  been  applied  to  finding  stress  intensity  fac- 
tors for  two  dimensional  problems  in  finite  regions,  including  applications 
of  the  method  of  Muskhelishvili , finite  el€anent  methods,  and  boundary  inte- 
gral techniques.  In  addition,  there  have  been  several  applications  of  the 

3 4 

crack  solutions  by  Williams  ’ as  the  elements  in  a series  expansion.  The 
first  element  in  this  series  contains  the  square-root  singularity  in  stress 
The  coefficient  of  this  term  is  defined  to  be  the  stress-intensity  factor. 
The  other  solutions  for  the  straight  crack  in  the  infinite  sheet  are  not 
singular,  but  provide  a denumerable  infinity  of  solutions  of  the  form 
a_(r,9)  = ^f^j(n,0),  n^2. 

A significant  advantage  in  using  the  Williams'  solutions  as  elements  in 
an  expansion  is  that  each  term  not  only  satisfies  the  field  equations,  but 
also  satisfies  the  boundary  conditions  at  points  near  to  the  singularity, 
i.e.,  on  the  edges  of  the  crack.  Consequently  the  errors,  whether  numer- 
ical or  due  to  truncation,  which  result  from  approximating  the  boundary 
conditions  on  a finite  domain  are  introduced  at  points  well  removed  from 
the  singularity. 

The  first  applications  of  the  series  of  Williams'  solutions  to  the 
finding  of  stress  intensity  factors  in  finite  sheets  appears  to  have  been 
made  by  Gross,  Srawley  and  Brown^’^’^,  in  which  the  necessary  coefficients 
were  found  by  boundary  collocation.  Stress  intensity  factors  for  edge 
cracked  rectangular  plates  in  tension  and  bending  and  for  other  geometries 

g 

were  obtained  in  this  manner.  More  recently,  Hulbert  and  Hopper  deter- 
mined expansion  coefficients  by  using  more  boundary  points  than  terms  in 
the  series  expansion,  and  determined  the  coefficients  so  as  to  minimize  the 
squared  error  on  the  boundary. 


1 


In  the  present  work,  the  Principle  of  Minimum  Potential  Energy  is  used 
to  develop  an  algorithm  for  finding  the  appropriate  coefficients  of  an  ex- 
pansion for  stresses  and  displacements  when  an  arbitrary  in-plane  traction 
is  applied  at  the  circumference  of  an  arbitrary,  edge-cracked,  finite 
domain.  It  is  also  shown  that  the  Principle  of  Minimum  Complementary 
Energy,  together  with  the  series  of  Williams'  solutions,  is  appropriate 
when  displacements  are  prescribed  on  the  edge  of  the  plate.  Finally, 
Reissner's  Principle  is  used  to  find  the  coefficients  of  the  expansion  when 
displacements  are  given  on  some  portion  of  the  boundary,  and  tractions  are 
prescribed  on  the  remainder. 


t 

i 

£ 

I 
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II.  Mef.hod  of  Solution 


Tlia  region  of  interest  i;;  depicted  p.repii icnlly  as  Figure  1.  A simply 
connected  elastic  sheet  of  negligible  thickness  (plane  stress)  is  assiaied 
to  contain  a straight  crack  of  length  a,  originating  at  an  edge.  For  con- 
venience, an  origin  of  coordinates  is  placed  at  the  crack  tip,  v;ith  crack 
boundaries  and  located  at  6 = + tt,  0^rj<a,  In-plana  tractions  are 
assuaed  to  be  prescribed  on  the  curve  S,  with  values  prescribed  to  be  zero 
on  and  S^. 
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It  has  been  shotai  by  Williams  ’ that  the  equations  of  elastostatics 
are  satisfied  for  the  infinite  domain  by  stress  functions  of  the  form 

X(r,0,X)  = r^'^^F(6;X)  (1) 

For  the  crack,  the  eigenvalues  X arc  X = ra/2,  for  all  integers  m.  For  the* 
coordinates  shown,  the  stress  function  may  be  divided,  for  convenience,  into 
an  even  and  an  odd  part.  The  even  stress  function  and  the  corresponding 
infinite  set  of  stresses  arc  given  in  Appendix  I.  For  each  n there  arc  two 
states  of  stress,  each  of  which  is  multiplied  by  an  arbitrarj'  coefficient 
and  each  of  which  satisfies  the  field  equations  and  the  traction  free  bound- 
ary condition  on  the  crack  edges.  For  simplicity  in  that  which  follows,  let 
us  designate  a superposition  of  M of  these  solutions  as: 

„ “V 
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where  dj^  = a^^,  d^  = -a^,  d^  = -a^,  d^^  - a^^,  etc. 

Our  purpose,  then,  is  to  develop  an  efficient 
the  set  of  M coefficients  which  give  tlie  best  appro 
xvithin  an  arbitrary  finite  domain,  using  tlie  first 


algorithm  for  findi 
xJmo.tiou  to  the  str 
M eigenfunctions. 


(2) 


'k; 

os.ses 

As 


noted  earlier,  both  collocation  and  the  Ee.thod  of  least  squares  have  been 
used  for  this  purpose.  Here,  wa  seek  to  use  an  ej-:trenial  principle  of  elasto- 


sta  tics . 


As  was  shotjn  by  Uillianis  ’ , displacements  corresponding  to  each  term 
in  the  stress  function  may  be  found  through 

"e  - -7?e'  + 37 

where  y is  the  shear  modulus,  v is  Poisson's  ratio,  o is  taken  to  be  v for 
plane  strain,  and  a = u/(l+v)  in  the  case  of  plane  stress.  The  function  i{/ 
which  corresponds  to  the  even  stress  function  and  the  resulting  displace- 
ments are  given  in  the  appendix. 

Thus,  for  any  m we  have  a set  of  displacements,  u^,  u^  which, 
upon  application  of  the  strain  displacement  equations  and  the  stress  strain 
laws,  w’ill  yield  the  stresses  *^r0^‘  the  class  oF  prob- 

lems  v;ith  prescribed  tractions  T on  the  boundary,  these  displacements  are 
kinetically  adnissable,  as  is  the  superposition. 


‘ M ^ 

j L * m 

2y  u^  m-1  jUg 

Such  displacements  are  therefore  appropriate  for  use  in  the  Principle  of 

9 

Minimum  Potential  Energy  , with  the  d to  be  chosen  to  achieve  the  minimu-m. 

m 


U(u^)dv  - 


f f * 

F.u.dv  - T.  u.dS 
11  J J r 1 


where 


U is  the  strain  energy  density 


F.  are  body  forces  and 
a 


(13) 


C = ICv  o + V )u  + (v  o + )u„  ids 

laq  J r rr  G 6r  r r rO  0 00  0 


= {T^.  (r,G)u^'^  Tg  (r,0)u^^}d£ 


The  desired  algorithm  for  choosing  the  coefficients  d has  now  been 
obtained.  For  any  preselected  number  of  terms,  M,  to  be  used  in  the  ex- 
pansion, the  elements  of  the  square  array  (which  can  be  shown  to  be 

symmetric)  are  evaluated.  The  prescribed  tractions  are  substituted  into 
Equation  14,  and  the  M vector  is  computed.  Post-multiplying  by  C ^ 

leads  to  the  coefficients  d . 

m 

For  purposes  of  determining  stress  Intensity  factors,  it  is  the  first 
of  these  coefficients  which  is  of  interest.  Let  v 

K = S,im  { /2Trr  a (9=0,r)  } ' (15) 

r->-0 

Only  the  first  term  survives  the  limit  process,  thus 

a^  = - (16) 

The  coefficients  d will  have  a linear  dependence  on  the  far  field 
m 

stress,  through  Equations  12  and  14,  thus  the  results  may  be  put  into  the 
usual  dimensionless  form  for  the  Mode  I stress  intensity  factor,  through 


a /i.a 
o 


wliere  a is  the  crack  length,  and  is  a suitable  measure  oF  the  far-field 
stress . 
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III.  Stress  Inteasity  Factors  for 
Rectangular  Sheets  in  Tension 


The  rectangular  sheet  xrith  edge  crack  loaded  in  uniaxial  tension  was 


selected  as  a suitable  problca  for  verifying  the  feasibility  of  the  method. 


The  problem  is  of  considerable  interest;  consequently,  previously  obtained 


results  are  available  for  comparison. 


We  consider  a rectangular  sheet  of  length  L and  width  b,  subjected  to 


a uniform  traction  T of  magnitude  c and  no  shear  on  the  ends,  the  other 

y O ’ 


faces  being  free.  The  geometry  of  the  problem  is  depicted  in  Figure  2. 


For  this  simply  connected  problem  x^ith  traction  boundary  conditions,  the 


exact  solution  should  be  independent  of  Poisson's  ratio.  However,  the 


present  method  introduces  Poisson's  ratio  tlirough  the  displacements.  A 


value  of  V = 0.3  was  used  in  the  calculations  here  reported. 


Since  each  of  the  stresses  and  o vanishes  on  6 = + tt,  inte- 

HH  rfl  — . ’ 


gration  over  that  portion  of  the  boundary  produces  no  contributions  to  the 


matrix  defined  by  F.quation  13.  For  the  problem  at  hand,  the  tractions 


and  C are  prescribed  to  be  zero,  thus  the  vector  D is  computed  by 


integration  over  only.  In  this  symmetric  problem,  only  the  upper  half 


plane  need  be  considered.  Since  the  boundary  tractions  on  the  plana  y = L/2 


are  prescribed  in  the  Cartesian  frame,  the  computation  of  the  invariant 


quantity 


“.-I 


T • u ‘ ds 


is  more  conveniently  performed  in  the  Cartesian  system,  i.e.. 


x=b-a 


D - T (u  ^ cos  C - u sin  0)  + T (u  ^ sin  0 + u„^  cos  0)dx(19) 
qjx  r 0 yr  0 


"k  k 

Here,  T =0  and  T = a . 

X y o 
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Since  evaluation  of  the  integrals  in  Equations  13  and  lA  in  closed 
fora  is  not  feasible  for  the  geometry  selected,  and  certainly  not  for  more 
complex  shapes,  a numerical  integration  using  S:impson's  rule  was  employed. 
Division  of  the  boundary  curve  into  200-'i00  points  v/as  found 

to  be  sufficient  to  keep  the  results  insensitive  to  the  number  of  points 
used.  The  approximations  introduced  through  those  numerical  integrations 
appeared  to  place  a limit  on  the  crack  depths  which  can  be  considered  by 
this  method.  Sjensitivity  to  the  number  of  integration  points  began  to 


appear  when  the  distance  from  the  crack  tip  to  the  nearest  boundary  was 
allowed  to  become  less  than  about  ten  times  the  distance  between  integra- 
tion points.  As  a practical  matter,  this  limitation  is  of  no  consequence, 
for  the  available  closed  form  solutions^^  for  very  shallow  and  very  deep 
cracks  in  strips  of  infinite  length  are  applicable  to  such  geometries. 

Selected  values  of  the  Mode  I Stress  Intensity  Factor  are  tabulated 
in  Table.  I for  various  crack  lengths  and  aspect  ratios.  A graphical  com- 
parison of  some  of  these  results  Is  given  as  Figure  3.  The  form  used  by 
Tada^^  was  found  to  be  particularly  convenient  for  this  comparison.  The 
use  of  a dimensionless  factor,  N,  defined  through 


K (a/b,L/b) 

N(a/b,L/b)  = — (1  - a/b)  ' 


o vTia 
o 


was  found  to  reduce  the  large  variations  evident  in  Figure  3 to  the  more 
readily  useable  results  presented  as  Figure  4. 

Numerical  values  obtained  at  small  a/b  (0.05  to  0.10)  appear  to  con- 
verge satisfactorily^^  to  the  knov.';>  value  for  a/b  o,  as  do  the  numerical 
values  obtained  at  large  a/b  (0.9  Lo  0.95).  V'aluer.  of  N obtained  for 
larg;>  I./b  (f./b  ~ 2,  3 arid  4)  give  satisfactory  agreement  v;ith  the  curve 


prusent'''d  by  Tnda  n.~.  a composite  of  previously  obtaiu.'d  v.ilues  ft»r  long 


7T  a - DIMENSIOMLESS  STRESS  INTENSITY 


0 .2  A .6  .8  1 


B/b  - CRACK  DEPTH 


Figurfi  3 


Stress  Intensity  Factors  for  Lc!ge  Cracked  Plates  of  Several 
Asuoct  Ratios 
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and  infinite  strips  by  Eowie  and  Keal  , Eiaery  and  Seyediu  , and  others. 

Expansions  using  from  eight  to  48  tenss  v/t-re  investigated.  In  most 
cases,  24  to  36  terms  were,  found  ade-iuate  to  give  convergence  to  three  sig- 
nificant digits  in  the  value  of  the  first  coefficient.  More  terns  ware 
generally  required  for  good  convergence  as  the  ratio  L/b  %vas  made  large 
and  as  the  crack  depth  approached  the  plate  edge.  From  results  such  as 
those  given  in  Table  II,  it  is  suggested  that  convergence  in  the  third 
significant  figure  \»'ill  only  be  achieved  if  the  most  rapidly  oscillating 
tern  has  half  period  less  than  the  smaller  of  and  (Figure  2),  i.e., 

M>j)  ^ 2ir.  Some  evidence  of  numerical  difficulty  began  to  appear  in  48  tern 
expansions.  These  were  attributed  to  the  large  variations  between  the 
magnitudes  of  the  elements  of  the  C matrix  which  occurs  in  such  cases. 

Thus,  the  method  is  not  well  suited  to  extremely  long  strips.  This  again 
is  of  no  consequence,  as  results  for  such  geometries  are  readily  obtained 
by  using  the  known  results  for  the  infinite  strip  and  St.  Venant's  Princip- 


The  determination  of  a single  set  of  coefficients  on  a CDC  6600  compu- 


ter with  a program  which  was  not  particularly  efficient  was  found  to  re- 

2 03 

quire  a run  time  of  approximately  0.3M  N x 10  sec,  tJhere  N is  the  num- 

P P 

ber  of  boundary  points  used  for  the  numerical  integration. 
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Table  II 


IV.  The  Influence  of  Notches 


The  npplicaf)  ility  of  the  procedure  here  presented  to  complex  boundar- 
ies was  used  to  investip.ate,  in  a limited  nannec,  the  influence  of  notch 
v/idth  on  the  stress  intensity  factor  for  cracks  originating  from  the  notch 
root.  The  geometry  considered  is  shown  in  Figure  5.  A triangular  notch 
of  depth  d and  width  2h  is  assumed  to  be  machined  into  the  edge  of  a rec- 
tangular tensile  sheet  of  width  b and  length  L.  A straight  crack  of  length 
a-d  is  assumed  to  originate  at  the  notch  root.  For  h = 0,  the  problem  re- 
duces to  that  of  an  edge  crack  of  depth  a.  The  numerical  integration  in 
Equation  13  was  accomplished  by  dividing  the  boundary  segment  into  two 
segments,  AB  and  BC.  Since  the  notch  faces  are  free  of  traction,  no  change 
in  the  vector  D results. 

q 

Numerical  results  were  obtained  for  a square  plate  (L/b  = 1)  with 
notch  depth  d/b  = 0.25.  ' The  distance,  a,  from  the  edge  of  the  plate  to  the. 
crack  tip  was  varied  from  a/b  = .3  to  a/b  = .5.  Dimensionless  stress  inten- 
sity factors  are  compared  in  Table  III  for  several  notch  widths.  It  is 
evident  from  the  results  that  notches  have  little  effect  on  the  stress  in- 
tensity factor  once  the  crack  is  well  established  (3  /d  0.2).  Thus,  it 
is  to  be  concluded  that  the  practice  employed  in  experimental  work  of  using 
a narrow  machined  cut  as  a crack  starter  and  then  considering  the  depth  of 
cut  to  be  a portion  of  the  crack  is  not  a source  of  signlficaut  error. 


V.  Mixed  Boundary  Value  Problems 
->-* 

If  displacements,  u , rather  than  stresses,  are  prescribed  on  the 
boundary  S2  (Figure  1)  and  the  tractions  are  prescribed  to  be  zero  on 
and  S^,  as  before,  a mixed  boundary  value  problem  results.  For  such  prob- 
lems, each  of  the  Williams'  solutions  is  a statically  admissable  system  of 
stresses  since  each  satisfies  equilibrium  and  the  traction  boundary  condi- 
tions. Thus,  any  superposition,  such  as  Equation  2,  is  also  statically  ad- 


missable and  may  be  used  in  the  Principle  of  Minimum  Complementary  Energy  . 


fff 

= W (a . . 

JJJ  c ij 


-V*  ->■ 

)dV  - u • T ds 


The  variation  of  it  is  taken,  with  the  stresses  treated  as  the  independent 

variables.  Variations  of  the  stresses, 

M 

6o.  . = y 6d  a. (22) 
X3  p^l  P 

satisfy  the  equilibrium  equations  and  the  traction  free  condition  on  S^+S^  ^ 
Thus,  setting  the  first  variation  of  the  Complementary  Energy  to  zero  leads 


j-  . M M 

{2pu.  - y d u.™}  * y 6d  v.a..^ds  = 0 

^ m=l  ^ p=l  P J 


The  same  choice  of  variation  coefficients  as  given  in  Equation  11  is 
apropos  ; thus  the  algorithm  for  determining  the  coefficients  of  the  ex- 
pansion is  seen  to  be 

M 

y d B = E ( 

^ , m mq  q 

m=l  ^ ^ 


where 


B “ [ V . ^ ds 

mq  j 1 J jx 


r.  * q 

E = 2-  u.  V.P.,  ds 

q 1 1 1^ 

^2 
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One  complication  arises  in  this  mixed  problem  which  is  not  present 
when  tractions  alone  are  prescribed.  Each  term  in  the  expansion  can  be  seen 
(Appendix  I)  to  have  zero  displacement  at  the  origin  of  coordinates.  Thus, 
the  crack  tip  is  automatically  specified  to  be  a fixed  point.  When  trac- 
tions are  prescribed  on  S this  is  of  no  consequence,  and  serves  to  remove 
the  rigid  body  displacements.  When  displacements  are  to  be  specified  on 
some  portion  of  S,  however,  the  crack  tip  must  be  left  free  to  translate. 
This  may  be  accomplished  by  including  the  two  possible  rigid  body  displace- 
ments as  additional  elements  of  the  expansion.  These  solutions,  in  polar 


coordinates,  are 


2n  u = d u = d cos  0 
r o r o 

2p  u„  = d u„°  = -d  sin  0 
9 o 0 o 


for  X axis  translation,  and 

2)j  u = d ,u  ° = d , sin  0 

' (28) 
2y  Ug  = = d_j^  cos  0 

for  y axis  translation.  Of  these,  only  the  former  (Equation  27)  is  symme- 
tric with  respect  to  the  x axis  and  appropriate  for  inclusion  with  the 
symmetric  solutions  given  in  the  appendix.  This  solution  represents  a 
rigid  body  motion:  thus  the  corresponding  stresses,  a o and  a 
all  vanish. 

A more  Interesting  class  of  problems  arises  if  the  boundary  S2  (Figure 
1)  is  allowed  to  be  divided  into  two  regions,  and  C^,  with  tractions 
prescribed  on  and  displacements  on  C^.  Neither  the  Principle  of  Minimum 
Potential  Energy  nor  the  Principle  of  Minimum  Complementary  Energy  can  be 
applied  to  such  problems,  for  the  Williams'  solutions  are  then  neither 
statically  admissable  nor  kinematically  admissable.  It  will  be  shown. 


however,  tliat  Rcissner's  Principle  can  be  used  to  develop  an  algorithm  for 
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detenaining  the  desired  set  of  coefficients. 


We  consider  the  functional 


J=  {W(c..)-F.u.}dV-  o..{e..-^(u.  ,+u.  .)}dV 

ij  11  J ±j  ij  2 i.j  j,i 

V V 

* f * 

- T.  • u.ds  - T.(u.  - u.  )ds 

J 1 1 J X 1 1 


where  S C and  S,^  = C..,  + S,  + S and  consider  stresses,  strains  and 
u u r T 1 3 

displacements  all  to  be  independent  variables.  We  let  the  displacements 
be  given  by  Equation  5,  the  stresses  by  Equations  2,  and  the  strains  by  an 
analagous  form.  Each  element  of  the  expansion  then  satisfies  the  equili- 
brium equations,  the  strain  displacement  equations,  and  Hookes'  law. 
Setting  the  first  variation  to  zero  leads  to 

6J  = 0 = [ T.  • 6u.ds  - [ T.*  • 6u.ds  - [ 6t.  ' (u.  - u *)ds 


0 = (T 


T.du.ds 

S 

u 

* f 

r - T ) • 6u.ds  - 6T.  ■ (u.  - u. 

i i 1 J 1 11 


Since  T = T * = 0 on  S.„  - = S,  + S_  for  each  element  of  the  expansion, 

i a T T 1 J 

Substituting  Equations  9 and  22,  and  making  the  same  choice  of  variation 
coefficients  as  specified  by  Equation  11,  we  find  that  the  desired  ex- 
pression for  the  expansion  coefficients  is 


r ^ 

I V . y d ."*u.*^ds 


- [ y d u.”''ds  =[  T.  u.*^ds  - [ 2p  u.  v.o..* 


VI.  The  Rectangular  Strip  with  Prescribed 
End  Displacements 


The  method  developed  in  the  preceeding  section  for  determining  co- 
efficients for  displacements  and  stresses  expressed  as  a series  of  Williams' 
solutions  can  be  applied  to  the  determination  of  stress  intensity  factors 

for  the  rectangular  sheet  with  two  free  edges,  prescribed  end  displacement, 

14 

and  an  edge  crack.  This  problem  has  been  considered  previously  by  Bowie  , 
who  used  a conformal  mapping,  and  by  Bradley  and  Kobayashi^^,  who  have  fit 
several  terms  from  the  Williams'  solutions  to  data  obtained  from  photoelastic 
experiments.  The  geometry  of  interest  is  shown  in  Figure  6.  Surfaces  S^, 

C , and  C are  free  of  traction  and  displacements  are  prescribed  on  C 

Ti  T2 

For  this  symmetric  geometry,  and  for  symmetric  loadings,  i.e.. 


u (x,L/2)  = - u (x,-L/2);  u (x,L/2)  = u (x,-L/2) 
y y X X 

only  the  upper  half  plane  need  be  considered.  Thus, 


mq 


I f I Q . / . m,  q,  , 

{v  a + v„a.  )u  ^ + (v  o . + v.a„„  )u„^}ds 

J r rr  0 0r  r r rO  0 00  0 


Ct  +C 
1 2 


(35) 


(36) 


F = - 2p 

q 


f * ^ GO 

{ (u  cos  0 + u sin  0)(v  0 + v„o„  ^) 

X y r rr  0 0r 


(37) 


+ (-U  sin  0 + u cos  0)  (v  a + v.o„„*^)}ds 
X y r r0  0 00 

and  the  d arc  obtained  from  Equation  33. 
m ^ 

Numerical  results  were  obtained  for  a uniform  extension  with  no  trans- 
verse displacements  on  the  ends  of  a sheet  of  length  T.  and  width  b,  having 
an  edge  crack  of  length  a.  The  geometry  is  shown  in  Figure  7.  The 
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prescribed  displacements  are: 

* 

u (x,L/2)  = u and 

y y (38) 

u^(x,L/2)  = 0 

Results  were  obtained  by  using  34  of  the  Williams'  solutions  and  th® 
rigid  body  mode  previously  discussed  as  the  elements  of  the  expansion.  The 
numerical  integration  was  found  to  require  300-400  points  on  the  boundary. 
Convergence  of  the  numerical  solution  was  found  to  be  similar  to  that  pre- 
viously discussed  for  prescribed  boundary  tractions. 

A convenient  means  of  reducing  the  Mode  I Stress  Intensity  Factor  to  a 


dimensionless  form  is  through  division  by  a^vna,  where 


2u 


is  the  uniform  traction  which  must  be  applied  on  y = + L/2  in  order  to 
produce  the  same  displacements  in  an  uncracked  sheet  of  the  same  material, 
subject  to  complete  transverse  constraint,  c =0.  Thus 

XX 


a Vva 
o 


(l-v2)  (-d  ) 

a 1 


Numerical  results  for  several  aspect  ratios  are  presented  as  the  dashed 
curves  of  Figure  8 for  the  case  v = 0.3. 

The  results  obtained  for  L/b  =1.0  appear  to  be  in  good  agreement  with 
those  plotted  by  Bradley  and  Kobayaski^^.  The  present  results  can  also  be 
compared  with  the  exact  result  for  a limiting  case.  For  an  clastic  strip 
of  infinite  width,  having  displacements  u^  on  clamped  edges  at  y = + L/2, 
Rice^^  has  shown  the  stress  Intensity  factor  to  be 


~ ^ ""y  L(l-v2) 


This  result  is  compared  with  the  present  predictions  in  Table  IV.  It 
should  be  noted  that  the  normalization  is  with  respect  to  square  root  of 
the  plate  width,  rather  than  the  square  root  of  the  crack  length,  in  order 
to  facilitate  the  comparison.  It  can  be  seen  that  the  stress  intensity 
factor  in  short  specimens  is  very  uniform,  and  close  to  that  given  by 
Equation  41. 

Symmetrically  located  center  cracks  in  strips  with  prescribed  end 

17  18 

displacements  and  free  edges  have  been  considered  by  Bowie  and  Isida 
Stress  intensity  factors  for  such  problems  were  also  found  to  be  nearly 
independent  of  crack  length  for  all  but  the  shortest  cracks.  Isida 's 
results  show  that  the  stress  intensity  factor  for  central  cracks  to  be 
qualitatively  the  value  given  by  Equation  41  for  L/b  < 1 and  a < c/2. 


L/b  = 0.5 

L/b  = 1.0 

L/b  = 2.0 

L/b  = 3.0 

a/b  =0.1 

0.251 

0.285 

0.295 

0.340 

a/b  = 0.2 

0.261 

0.352 

0.440 

0.473 

a/b  = 0.3 

0.261 

0.371 

0.508 

0.573 

a/b  = 0.4 

0.266 

0.375 

0.534 

0.639 

a/b  = 0.5 

0.264 

0.375 

0.539 

0.670 

a/b  = 0.6 

0.264 

0.375 

0.533 

0.673 

a/b  = 0.7 

0.266 

0.373 

0.523 

0.653 

a/b  = 0.8 

0.265 

0.370 

0.507 

0.619 

a/b  = 0.9 

0.260 

0.370 

0.510 

0.606 

b = <x> 

0.269 

0.381 

0.538 

0.659 

exact 
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VII.  Calibration  of  the  Fatigue  Crack  Life  Gauge 


19  20  21 

It  has  been  proposed  ’ ’ that  a thin  sheet  of  metal  containing  a 

crack  (the  gauge)  be  attached  to  a structure  and  the  growth  of  the  crack 

in  service  monitored  as  a means  of  predicting  the  rate  of  crack  growth  in 

19  20  21 

the  structure.  Both  edge  cracks  , and  central  cracks  ’ have  been 
suggested,  although  the  edge  crack  would  appear  to  offer  advantages  of 
ease  in  fabrication.  It  may  also  be  noted  that  the  observed  growth  rate 
of  the  gauge  crack  provides  a much  needed  means  of  quantifying  the  dam- 
aging aspects  of  the  service  load  spectrum. 

The  gauge  and  plate  are  shown  in  Figure  9.  A perfect  bond  is  assumed 
on  line  AB  and  CD,  with  gauge  and  plate  not  interacting  elsewhere.  A 
straight  edge  crack  of  length  a is  assumed  to  be  located  on  the  axis  of 
symmetry  of  the  gauge.  It  is  assumed  that  the  far-field  stresses  in  the 
plate  are  a uniaxial  tension,  a^.  The  state  of  strain  in  the  plate  is 
uniform  if  forces  generated  in  the  gauge  produce  negligible  deformations 

in  the  plate,  i.e.,  if  (E  t /E  t ) <<  1.  In  this  case,  the  displacement  of 

g 8 P P 

AB  with  respect  to  0 is 

Uy(x,L/2)  = Uy*  = (|)  (42) 

Points  on  the  line  AB  also  undergo  displacements 

u^(x,L/2)  = u/  = - VpX  (43) 

due  to  the  Poisson  contraction  of  the  plate  to  which  the  gauge  is  presumed 

to  be  perfectly  bonded.  Here  E and  v are  Young’s  modulus  and  Poisson’s 

P P 

ratio  forthe  plate,  respectively. 

These  displacements,  together  with  the  condition  that  the  edges  C„ 

and  C are  free  of  traction,  furnish  the  boundary  conditions  for  the 
^2 

fatigue  crack  life  gauge.  Substituting  the  displacements  into  Equation  37 


i 

1 


I 


permits  the  coefficients  of  the  expansion  for  the  stresses  and  displace- 
ments in  the  gauge  to  be  determined. 

Thirty-five  elements  were  again  used  in  the  expansion,  with  the  numer- 
ical integration  of  Equation  32  being  carried  out  at  300-400  boundary  points. 
The  mode  I stress  intensity  factor  is  conveniently  made  dimensionless 
through  division  by  where  is  the  far  field  stress  intensity  in 

the  plate.  We  find 


K, 


2uy*  _ 

E (-J ) Aa 

p L 


(44) 


Dimensionless  stress  intensity  factors  for  a gauge  mounted  on  a plate 
are  given  in  Figure  8,  where  they  may  be  compared  with  those  previously  given 
for  the  case  of  complete  transverse  constraints.  A Poisson's  ratio  of  0.3 
was  used  in  these  computations.  It  is  of  interest  to  note  that  the  differ- 
ences are  not  significant,  generally  being  within  the  1-v^  stiffening  effect 
due  to  transverse  constraints.  The  Mode  I stress  intensity  factor  is  essen- 
tially determined  by  the  displacements  perpendicular  to  the  plane  of  the 
crack,  and  is  not  sensitive  to  displacements  parallel  to  the  crack  plane. 

The  rigid  body  motion,  allowed  through  Inclusion  of  a rigid  iiorizontal 
translation  as  one  of  the  solutions,  was  found  to  have  a significant  ampli- 
tude. The  displacement  of  the  crack  tip  in  the  direction  of  the  crack  can 
be  seen  from  Figure  10  to  be  a significant  fraction  of  the  end  displacement 
(longitudinal).  These  displacements  may  be  sufficient  to  lead  to  difficul- 
ties in  attempts  to  view  (photographically)  the  crack  tip  while  a cyclic 
test  is  in  progress. 

Stress  intensity  factors  for  the  fatigue  crack  life  gauge  with  0.1  ^ 
a/b  ^ 0.9  are  presented  in  Figure  11  after  a different  normalization.  It 


32 


X 

! 


r 


can  be  seen  that  the  stress  intensity  factor  in  gauges  of  low  aspect  ratio 

does  not  vary  significantly  as  the  crack  length  changes.  Thus,  a gauge 

crack  under  constant  load  will  grow  at  a relatively  constant  rate,  and  the 

rate  of  change  of  crack  length  wil]  provide  a measure  of  instantaneous  load 

spectrum  severity  which  is  essentially  independent  of  instantaneous  crack 
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length.  A constant  growth  rate  under  constant  load  has  been  observed 

Verification  of  the  results  may  be  obtained  through  consideration  of 
two  special  cases.  Rice^^  has  given  a stress  concentration  factor  for  the 
centrally  cracked  infinite  strip  with 

* 

u (x,  + L/2)  = + u on  y = + L/2 

(45) 

T (x,  + L/2)  = 0 
xy  — 

as  Kj  = a</hJ2  where  (46) 

a = E u *2/L  (47) 

g y 

The  boundary  conditions  for  the  present  problem  (Equations  42  and  43)  lead 
to  the  same  energy  density  at  large  distances  in  front  of  the  crack  as  do 
those  of  Equation  45.  Consequently,  the  argument  given  by  Rice  may  be  used 
to  deduce  that  the  stress  intensity  must  be  the  same.  After  some  manipula- 
tion, we  find: 

(— ) = /—  (48) 

Eg  n/ 

It  may  be  seen  from  Figure  11  that  this  prediction  (arrows)  gives  surprisingly 
good  agreement  for  L/2b  < 1,  even  though  Equation  48  is  actually  appropriate 
only  for  lim.  {L/2b  »•  0}. 

A second  limiting  case  is  appropriate  for  a/b  ->■  0 and  a/L  ->■  0.  In 
this  case,  the  stress  intensity  factor  should  be  that  of  an  edge  cracked 


infinite  sheet,  or 


(/)  = 1.12 

P 
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The  present  results  show  good  agreement  with  this  prediction  (dashed  curve. 
Figure  11)  for  small  a/b  and  long  strips. 

A properly  designed  gauge  should  display  neither  extremely  high,  nor 
extremely  low,  sensitivity.  That  is,  a gauge  length  which  is  entirely 
consumed  during  insignificant  crack  growth  in  the  structure  is  a poorly  de- 
signed gauge,  as  is  a gauge  which  shows  little  growth  during  significant 
(and  dangerous)  growth  in  the  structure. 

At  some  critical  point  in  the  structure  a crack  of  length  a^  may  be 
presumed  to  be  growing  in  accord  with  the  Paris  Law,  with  growth  rate  dom- 
nated  by 

K = a Aa  M (50) 

Ig  P s 

where  M is  typically  greater  than  1.  In  the  gage,  the  crack  growth  rate 
is  essentially  that  given  by  Equation  45,  or,  for  very  short  cracks,  less. 
Hence, 

and  we  have 


A 

where  M < 1 and  decreases  as  the  cracks  grow.  Thus  the  initial  ratio  of 

gauge  crack  growth  rate  to  structure  crack  growth  rate  can  be  specified 

through  the  selection  of  guage  length  and  the  modulus  of  the  gauge.  An  ini- 
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tial  crack  length  greater  than  L/2  has  been  suggested  as  being  necessary 

to  assure  a constant  growth  rate.  The  results  given  as  Figure  11  suggest 

that  this  criterion  may  be  unnecessarily  conservative,  and  that  a ratio 

a /L  > 0.25  would  suffice, 
o 

As  both  cracks  grow,  the  ratio  of  stress  intensity  factor  in  the  gauge 


K - r -S-  /I 

p E J 2 
g P 


/ L , * 


C £ 

K~  = E /2na 
I p V s 

s 


M 


(51) 


(52) 


to  that  in  the  structure  will  diminish,  the  crack  in  the  structure  will 
begin  to  grow  more  rapidly,  and  the  gauge  will  loose  it's  sensitivity. 


It  may  therefore  be  necessary  to  use  more  than  one  crack  life  gauge  on  a 
structure,  with  one  being  a durability  gauge  designed  for  optimum  sensitiv 
ity  as  cracks  in  the  structure  approach  the  length  at  which  crack  main- 
tenance is  required,  and  the  other  gauge,  a safety  gauge,  being  designed 
for  optimum  sensitivity  as  cracks  in  the  structure  approach  the  critical 
crack  length. 
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VIII.  DISCUSSION  AND  CONCLUSIONS 


f 


i 


An  eigenfunction  expansion  was  used  to  describe  the  stresses  and  dis- 
placements in  a finite  elastic  sheet  containing  a single  straight  crack 
which  intersects  an  edge.  The  set  of  solutions  obtained  by  Williams  for 
the  infinite  sheet  with  crack  were  used  as  elements  of  the  expansion,  with 
coefficients  chosen  by  means  of  the  energy  principles  of  elastostatics . 

The  procedure  generates  an  approximate  solution  which  satisfies  the  field 
equations  everywhere  and  satisfies  the  traction  free  condition  on  the  crack 
face.  l-Zhen  tractions  or  displacements  are  prescribed,  a stationary  value 
of  potential  energy  or  complementai^r  energy  is  achieved.  Such  solutions 
approximate  the  boundary  conditions  elsewhere  as  well  as  can  be  accomplished 
with  a given  number  of  terms.  The  use  of  Reissner ’s  Principle  was  found 
to  make  possible  the  construction  of  solutions  when  mixed  boundary  condi- 
tions are  given. 

Although  the  numerical  results  do  not  suggest  any  difficulties  with 
convergence,  the  set  of  functions  employed  has  not  been  demonstrated  to  be 
complete.  No  multiply-connected  domains  were  considered.  The  boundary  may 
have  arbitrary  shape,  although  it  is  to  be  expected  that  boundaries  with 
re-entrant  corners,  other  than  the  single  crack,  might  lead  to  difficulties. 
No  attempt  was  made  to  investigate  the  convergence  at  singularities  intro- 
duced by  abrupt  changes  in  boundary  conditions. 

The  functions  are  not  orthogonal,  in  that  forces  of  the  pth  mode  do 
work  in  acting  through  the  displacements  of  the  qth  mode.  Consequently, 
the  array,  C,  is  not  diagonal  and  the  addition  of  elements  to  the  series 
leads  to  changes  in  all  coefficients.  Convergence  was  assumed  to  exist 
when  the  addition  of  terms  to  the  expansion  did  not  change  the  coefficient 


of  the  first  term.  The  most  rapid  convergence  was  found  to  occur  when  the 
boundary  points  at  which  the  prescribed  tractions  were  approximated  were 
all  at  about  the  same  distance  from  the  crack  tip.  Conversely,  more  terms 
were  found  to  be  required  for  crack  tips  near  boundaries,  and  for  large 
aspect  ratios.  A limitation  in  the  number  of  terms  which  can  be  used, 
other  than  one  caused  by  computer  storage  capacity,  was  encountered  and  may 
limit  the  degree  of  accuracy  achievable  with  the  method.  Eigenfunctions 
having  rapid  angular  oscillations  are  required  for  some  geometries,  and  the 
large  powers  of  radius  which  occur  in  such  cases  lead  to  large  variations 
in  the  magnitudes  of  the  elements  of  C.  This  leads  to  truncation  errors 
in  the  solution  of  the  system  of  equations. 

For  the  rectangular  sheets  considered  as  examples,  these  limitations 
precluded  obtaining  accuracy  results  for  the  stress  intensity  at  crack  tips 
near  the  boundary,  or  for  long  strips.  Stress  intensity  factors  for  the 
former  case  can  readily  be  obtained  from  the  well  known  results  for  infi- 
nite sheets.  Results  for  long  strips  with  prescribed  tractions  are  not  re- 
quired, as  it  was  found  that  a strip  with  an  aspect  ratio  of  three  is,  for 
all  practical  purposes,  infinite  in  length.  Long  strips  with  prescribed 
end  displacements  are  somewhat  more  complex,  although  a procedure  for 
estimating  the  stress  intensity  factor  from  the  stress  intensity  factor 
and  compliance  of  short  strips  is  being  developed. 

Results  obtained  for  the  edge  cracked  rectangular  sheet  have  several 


Important  Implications  for  the  design  of  fatigue  crack  life  gauges.  These 
are:  (1)  The  edge  cracked  configuration  appears  more  desirable  than  the 

center  cracked  configuration,  for  a shorter  initial  crack  length  is 
necessary  in  order  to  achieve  a constant  stress  Intensity.  This  provides 
more  than  twice  the  useful  crack  growth  region.  (2)  Nearly  uniform  stress 


intensities  are  achievable  in  square  sheets.  (3)  It  may  prove  necessary 
to  use  gauges  which  are  of  a different  material  than  the  structure  in  order 
to  tailor  the  relative  growth  rate.  Materials  with  different  moduli,  as 
well  as  crack  growth  parameters,  should  be  considered.  (4)  More  than  one 
crack  gauge  may  be  required  on  any  structure,  each  having  a different 
sensitivity,  as  it  may  not  be  possible  to  design  a gauge  having  the  desired 
sensitivity  over  the  entire  range  of  structural  crack  growth. 

Although  no  anti-.:ymmetric  problems  were  considered,  it  is  to  be  ex- 
pected that  the  formal  computation  of  Mode  II  stress  intensity  factors  for 
straight  cracks  in  finite  domains  can  be  performed  in  the  same  manner.  The 
even  functions  given  in  the  appendix  must  then  be  replaced  with  the  appro- 
priate odd  functions,  and  suitable  traction  applied. 

In  addition  to  its  utility  in  the  determination  of  stress  intensity 
factors,  the  method  can  readily  be  applied  to  the  determination  of  stresses, 
strains  and  displacements  in  two-dimensional  regions  having  re-entrant  cor- 
ners. For  such  problems,  the  appropriate  eigen-functions  are  those  corres- 

3 

ponding  to  the  eigen-values  given  by  Williams  for  infinite  wedges  of  ver- 
tex angle  greater  than  180°. 
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Append L>:  I 

The  syniinetric  state  of  stress,  i.e.,  and  even  in  6,  odd. 


rr 


re 


is  generated  by  the  even  stress  function 


Xg  = (“1)  ^2n-l^  ^ " cos(n-  -)6  + cos(n+’i)0} 

+ (-l)^a2^r^^^  { - cos(n-l)0  + cos(n-M)O) 


(1) 


The  corresponding  stresses  are  found  frora  the  usual  relationships 


o = 
rr 


= 


_1  1 ^ 
r^  302  r 3r 


00  0r 


r0 


r 9r39  r^  38 


(2) 


(3) 


(4) 


= (-I)*'  ^r’^  ^^~a2^_j^(n-3/2)  (n-1/2)  {sin(n+ij)0  - sin(n-3/2)0} 
+ (-l)^r*^  ^a^^Cn)  { (n+l)sin(n+l)9  - (n-l)sin(n-l)9} 


(5) 


, 1^n-l  n-3l2 

°rr  - <-« 


, |'t(n-V2)^  - (i>2Ji)]cos(n-3/2)8\ 


I- 


(p.-3/2)  (n-’-2) cos  (n+^i)  0 


I / -1  xU  n 1 
+ (-1)  r a 


2n 


• I [(n-1)  - (n+1) ] cos (n-1) 0 

/ 

- (n+L)n  cos(nM)0 


(b) 


o„^  = (-.1)"  ^r"  (n-'-j)  f cos(n+’;)8-  cos(n-3/2)0  } 


09 


+ ^aj,^(n)(nH)  {cos(ii+l)fi  - cos(u-l)9} 


(7) 


It  i.s  easily  verified  that  the  stresses  a and  o vanish  on  0 = + v for 

0 J rO 


all  n. 


For  each  ij)^,  we  nay  compute  the  harnonic  function  required  to  gener- 


ate the  displacements  by  solving  tha  differential  equation 


o ^ ^ G 

^ aF-  > 


(8) 


The  result  is: 

0 = (-1) 
c 

+ (-1) 


n-3/2 

n-1 

r 

n~h  ^2n-l 

n r 

n 

2n  L 

• r(n-3/2)^-  - (nj^ 


n-3/2 


sin(n-3/2)0 


2 2 
- (n+1)^ 


siii(n 


-i)ej 


The  corresponding  displacements,  using  Equations  3 and  4 are 

, ( [(n+Jg)  + -^^{(n-3/2)^  - (n-?2)^}]cos(n-3/2)0 

2n-l  ! 

^ - (n-3/2) cos (n+4)0 


2y  = ( 


I t 1 \ ^ ^ 

+ (-1)  r a 


[(n+l)  + — ■ {(n-1)^  - (n+l)^}]cos(n-l)9 

n 


2n 


- (n+l) cos (n+l) 9 


and 


„ , ,.n-l  n~h 

2\i  u„  = (-1)  r a 


[-(n-3/2)  + {(n-3/2)^  - (n+?g)^}]sin(n-3/2)9' 


2n-l 


+ (n-3/2) sin (n+%) 9 


(11) 


+ (-l)\"a 


[-(n-1)  + {(n-1)^  - (n+1)^}] sin(n-l)0 


2n 


+ (n+l) sin (n+l) 9 
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